Particle production in p — p collisions and prediction for LHC energy 
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We analyze recent data on particle production yields obtained in p — p collisions at SPS and RHIC 
energies within the statistical model. We apply the model formulated in the canonical ensemble and 
focus on strange particle production. We introduce different methods to account for strangeness 
suppression effects and discuss their phenomenological verification. We show that at RHIC the 
midrapidity data on strange and multistrange particle multiplicity can be successfully described by 
the canonical statistical model with and without an extra suppression effects. On the other hand, 
SPS data integrated over the full phase-space require an additional strangeness suppression factor 
that is beyond the conventional canonical model. This factor is quantified by the strangeness satu- 
ration parameter or strangeness correlation volume. Extrapolating all relevant thermal parameters 
from SPS and RHIC to LHC energy we present predictions of the statistical model for particle yields 
in p — p collisions at y/s = 14 TeV. We discuss the role and the influence of a strangeness correlation 
volume on particle production in p — p collisions at LHC. 
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I. INTRODUCTION 

Results from the Large Hadron Collider (LHC) have a 
potential to usher in a new era of discoveries and insights 
into particle physics. In order to fully appreciate these it 
is important to have a clear understanding as to what we 
expect will happen at the LHC without new physics and 
new phenomena involved. This requires an extrapolation 
of trusted models to LHC energies. In this paper we use 
the statistical model for particle production which has 
been around since more than half a century [3, 0, Q and 
discuss its predictions for LHC energy. 

The statistical model has been very successful in de- 
scribing hadron yields in central heavy-ion collisions [U, 
HI H 0] • It has provided a useful framework for describing 
ccntrality and system size dependence of particle produc- 
tion in heavy-ion collisions [!,[§]. It has also led to new 
insights about its applicability in small systems like p — p 
and even e + - 



scattering 0, E3, El 0, ■ 

In the following we concentrate on particle production 
in elementary collisions and discuss the predictions of the 
statistical model for particle yields in p — p collisions at 
the LHC at *Js = 14 TeV center of mass energy. 

Particle production calculated within the statistical 
model is quantified by a set of thermal parameters, these 
are: the temperature, the volume and the set of chem- 
ical potentials which arc related to conserved charges. 
To make any predictions for particle production at LHC 
energies one needs methods to extrapolate the thermal 
parameters to higher energies. In this paper we first an- 
alyze the recent experimental data obtained mp—p colli- 
sions at SPS and RHIC and then we extrapolate these to 
LHC beam energies. This allows us to make predictions 



for different particle yield ratios at the LHC as well as 
discuss their dependence on the values of extrapolated 
parameters. A very brief account of some of our results 
was given in In this paper we present a complete 

description of the extrapolation method which is based 
on our new analysis of the p — p data at SPS and RHIC 
energies. 

The paper is organized as follows: In Section |TT] we 
summarize the main features of the canonical statistical 
model. In Section Hill we present the analysis of SPS and 
RHIC data obtained in p — p collisions. In Section llVl we 
introduce the extrapolation of the model to the LHC en- 
ergy and discuss its predictions for particle production in 
p—p collisions. In the final section we present conclusions 
and summarize our results. 



II. STATISTICAL MODEL 

Systematic studies of particle yields extending over 
more than a decade, using experimental results at dif- 
ferent beam energies, have revealed a clear underlying 
freeze-out pattern for particle yields in heavy-ion colli- 
sions A detailed comparison of different freeze- 
out criteria was made in [16T ] and we followed the one 
used previously [ItJ to extrapolate the temperature T 
and baryon chemical potential [Ib to the LHC energy. 
The expected hadron multiplicity ratios can be calcu- 
lated directly from this extrapolation if the grand canon- 
ical (GC) description of charge conservation is adequate. 
This is because, in the GC ensemble, any particle multi- 
plicity ratio is uniquely determined by the values of the 
temperature and the chemical potentials. 
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FIG. 1: Particle ratios as a function of the volume radius 
R. The temperature T = 170 MeV and the baryon chemical 
potential [ib = 1 MeV were chosen according to the thermal 
conditions expected to be valid at LHC energy. All ratios are 
normalized to their grand-canonical values. 
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A. Canonical suppression 



FIG. 2: Integrated particle yields from p-p collisions at 
a/I = 17.3 GeV |H HI, [13, [H| together with different model 
predictions: (a), (b) and (c) introduced in the text. The lower 
panel shows the deviation of the model fits to data. 



The usual form of the statistical model in the grand 
canonical ensemble formalism cannot be used when ei- 
ther the temperature or the volume or both are small, as 
a rule of thumb one needs VT 3 > 1 for a grand canonical 
description to hold [H, [lj|. Furthermore, even if this 
condition is met, if the abundance of a subset of parti- 
cles carrying a conserved charge is small, the canonical 
suppression still appears even though the grand canon- 
ical description is valid for the bulk of the produced 
hadrons. There is by now a vast literature on the subject 
of canonical suppression and we refer to several articles 
(see e.g. @, Mm)- 

The effects of canonical suppression arc illustrated in 
Fig. [1] which shows particle ratios of strange and multi- 
strange hadrons to pions normalized to the values in 
the grand canonical limit as a function of the radius 
of the system. The smaller the volume and the larger 
strangeness content of the particle the stronger the sup- 
pression of the yield of strange particles . This has been 
discussed in great detail in [ljj, [22| . 

The analysis of the variations of particle ratios with 
the size of the system, e.g . via the number of partici- 
pants, at SPS revealed [13| that the experimental data 
show stronger suppression of strange-particle yields than 
that expected in the canonical model [HI, [24[ ■ Con- 
sequently, an additional suppression effect had to be in- 
cluded in order to quantify the observed yields of strange 
particles. 



B. Strangeness correlation volume 



One possible explanation for the failure of the canon- 
ical corrections is that strangeness can be conserved ex- 
actly in a small subvolume of the fireball, thus leading 
to a stronger canonical suppression, as seen in Fig. [TJ 
even though the strange particles within such subvol- 
ume are taken as being in chemical equilibrium. A 
modification of the statistical model was formulated 
in [l3l [2^ that allows to quantify this extra suppres- 
sion by the strangeness correlation volume (cluster size) 
within which the strangeness is conserved exactly. An 
alternative method was proposed in [Io| to include in the 
canonical statistical model an additional factor 75 that 
accounts for possible deviations of strange particle abun- 
dance from their chemical equilibrium distribution value. 

In the next section we discuss how to choose the size 
of the cluster or 75 in the canonical model to describe 
experimental data on strange particle production in p — p 
collisions. First, we use data from SPS and from RHIC 
energies and next we discuss the possible extrapolation 
of the model parameters to make predictions for particle 
production in p — p collisions at LHC energies. 
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III. DATA ANALYSIS 

At the top SPS energy data are available for 
7T, K, Kg p and A yields integrated over the full phase- 
space [25], HH, H3, HI] in p — p collisions. In Fig. [2] we 
show the comparison of the statistical model with these 
data for three different implementations of the canonical 
statistical model labelled by (a), (b) and (c). In all cases 
the conservation of the electric charge and baryon num- 
ber is formulated in the GC ensemble and is thus con- 
trolled by the corresponding chemical potentials. How- 
ever, strangeness conservation is always treated canoni- 
cally. We consider three different models: 

(a) the strangeness suppression, at fixed T and //_b, is 
controlled by the volume parameter which coincides 
with the system size, 

(b) an additional strangeness suppression is introduced 
through the 7s-fugacity factor, 

(c) the canonical suppression is controlled by the clus- 
ter volume which can be smaller than the system 
size. 

From the comparisons shown in Fig. [5] it is clear the 
model (a) fails to describe the p — p data. Although, 
the canonical description of strangeness production is in- 
cluded in this model, the strange particle ratios exhibit 
large deviations from data. The discrepancies seen in 
Fig. [2] for model (a) could be even larger for multistrangc 
particles. The reason of these discrepancies is due to the 
fact, that the volume of the system, fixed from the pion 
yields, is already too large to imply strangeness suppres- 
sion. From the \ 2 fit to SPS data one gets, in model (a) 
the system size radius R ~ 1.3. Consequently, as seen 
in Fig. at this value of R the canonical suppression 
of single-strange particles is indeed negligible. Thus, for 
S = ±1 particles the model (a) is almost equivalent to 
the GC treatment of strangeness conservation. 

From the above discussion it is clear, that at the SPS 
the strangeness suppression due to canonical effects alone 
is not sufficient to describe the An data. In models (b) 
and (c) we have included an additional suppression of 
strange particle phase-space by introducing either the 75 
factor or the strangeness correlation volume. From the 
comparisons of the model with data shown in Fig. [2] one 
sees that both these models describe the SPS data quite 
well with similar values for T and /j,b (see Table [J). The 
particle yields calculated in these models are summarized 
in Table HT1 The only essential difference between model 
(b) and (c) appears in the yield of <p meson which in both 
cases is inconsistent with data as seen in Table [TT1 [33j . 
for a detailed discussion see Ref. fl3j ] . 

In order to quantify the change of thermal parameters 
with collision energy we also analyze data on particle pro- 
duction in p — p collisions at RHIC. Here, particle yields 
are only available around midrapidity [30. l3l|. However, 
contrary to SPS the yields of multistrangc particles are 



TABLE I: Statistical model parameters, extracted from the 
comparison of model (c) (upper table) and model (b) (lower 
table) with experimental data on p-p collisions. Full phase- 
space (4-7T integrated) data at ^/s = 17.3 GeV and midrapidity 
densities at yfs = 200 GeV were used. 



v^(GeV) 


T(MeV) 


MMeV) 


i?(fm) 


Rc(fm) 


17.3 


169L16 


225L45 


1.4L0.3 


1.0L0.3 


200 


165L4 


14L18 


1.2L0.1 


1.3L0.1 


v^(GeV) 


T(MeV) 


MMeV) 


R(im) 


7s 


17.3 


164±9 


211L27 


1.5±0.2 


0.7L0.1 


200 


165L4 


14L18 


1.2L0.1 


1.1±0.1 



TABLE II: Particle yields (4-7T integrated) in minimum bias p- 
p collisions at <Js = 17.3 GeV from Ref. gl, [H [H, @ and fit 
results from models (b) and (c). The cj> yield (below the line, 
Ref. [2!|) was not included in the fit. The values given here 
are model predictions using the thermal parameters extracted 
from the analysis of the remaining data. 



Particle 


Data 


Fit (b) 


Fit (c) 




3.02 ± 0.15 


3.01 


3.01 


7T+ 


2.36 ± 0.11 


2.32 


2.32 


K~ 


0.258 ± 0.055 


0.268 


0.270 


K+ 


0.160 ± 0.050 


0.162 


0.166 


A 


0.116 ± 0.011 


0.119 


0.120 


A 


0.0137 ± 0.0007 


0.0135 


0.0134 


K° s 


0.18 ± 0.04 


0.21 


0.20 


P 


0.0400 ± 0.0068 


0.0464 


0.0469 





0.0120 ± 0.0015 


0.0271 


0.0598 



also available allowing for more complete verification of 
thermal models used in the analysis at SPS. In Fig. [3] we 
compare the predictions of models (a), (b) and (c) with 
RHIC data. The particle yields used in this study and 
the yields calculated in the model are summarized in Ta- 
ble [TH] The resulting thermal parameters for models (b) 
and (c) are indicated in Table HI The STAR Collabora- 
tion published the statistical model analysis of p — p data 
restricted to yields of pions, kaons and (anti)protons [30| . 
using the undersaturation factor 75. Our analysis agrees 
with Ref. if one uses the same data set. With the hy- 
pcrons included in the analysis there is a slight increase 
in 7s whereas the temperature and chemical potential 
stay almost the same. 

From Fig. [3] and Table IIIII one sees that all models 
describe RHIC data within two standard deviations. At 
RHIC the pion yield at midrapidity is by more than a fac- 
tor of two lower than at SPS in 4-7T, resulting in the corre- 
sponding decrease of the volume parameter and stronger 
strangeness suppression. In the context of the consid- 
ered RHIC data it is rather difficult to definitely verify, 
that there is a need for an extra strangeness suppression 
effects going beyond the one already included in the stan- 
dard canonical model. This is clear from Table HTT1 where 
in model (b) the 75 ~ 1 and in model (c) the cluster and 
the system volume coincides within errors. 
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FIG. 3: Midrapidity particle densities from p-p collisions at 
^/s = 200 GeV from STAR [13, HH together with the model 
fits as in Fig. [2] The lower panel shows the ^-deviations of 
the model fits to data. 



TABLE III: Particle midrapidity densities in minimum bias p- 
p collisions at y/s — 200 GeV from Ref . (3(| [3l| and fit results 
from models (b) and (c). Charged kaons were analysed with 
two different techniques and both measurements are included 
in the fits. The cj> yield (below the line, Ref. [13]) was not in- 
cluded in the fit. The values given here are model predictions 
using the thermal parameters extracted from the analysis of 
the remaining data. 



Particle 


Data 


Fit (b) 


Fit (c) 




1.42 ± 0.11 


1.28 


1.26 


7V + 


1.44 ± 0.11 


1.31 


1.30 


K~ 


0.145 ± 0.013 


0.148 


0.148 


K+ 


0.150 ± 0.014 


0.153 


0.153 


V 


0.113 ± 0.010 


0.125 


0.124 


V 


0.138 ± 0.013 


0.153 


0.151 


A 


0.0351 ± 0.0039 


0.0272 


0.0279 


A 


0.0385 ± 0.0042 


0.0317 


0.0324 




0.0029 ± 0.0013 


0.0016 


0.0017 




0.0026 ± 0.0011 


0.0018 


0.0019 


fi" + Q + 


0.00034 ± 0.00026 


0.00017 


0.00018 


K~ (kink) 


0.137 ± 0.013 


0.148 


0.148 


K+ (kink) 


0.140 ± 0.014 
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K% 


0.134 ± 0.014 
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0.0180 ± 0.0037 


0.0391 


0.0317 
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FIG. 4: Cluster radius Rc as a function of energy, extracted 
from the present analysis of SPS data (triangle) and RHIC 
data (star). Also shown is the previous analysis [l3|] of a 
smaller set of SPS data (square). The lines illustrate possible 
evolutions towards LHC energies as discussed in the text. 



In the following, we concentrate on the statistical 
model predictions for particle yields at LHC energy. In 
view of the results obtained at SPS and at RHIC, we limit 
our attention only to model (c) and its extrapolation to 
LHC energies. 



IV. PARTICLE RATIOS IN p 
AT LHC 



p COLLISIONS 



The extrapolation of particle ratios to LHC energy re- 
quires estimates of the temperature, the chemical poten- 
tial and the cluster volume. From our analysis made at 
SPS and RHIC energies (summarized in Tables [H] and 
IIIip it is clear that no variation in the temperature is ex- 
pected between SPS, RHIC and LHC. A strong decrease 
of [Lb from SPS to RHIC seen in our results, together 
with the previous systematics on the beam energy de- 
pendence of \ib obtained from freezeout conditions in 
heavy- ion collisions indicate that all chemical po- 
tentials should be very small at LHC. In the following, 
we use T ~ 170 McV and [i B ~ 1 MeV from Ref. [Hi as 
appropriate thermal parameters at LHC. This tempera- 
ture can be interpreted as the critical temperature in the 
QCD phase transition and/or as the Hagedorn limiting 
temperature of the hadronic medium. 

The only parameter that remains largely uncertain in 
extrapolating from SPS and RHIC to LHC is the size 
of the cluster described by the radius Rc quantifying 
the strangeness suppression. Two limiting cases for the 
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FIG. 5: Plot Predictions for various particle ratios using dif- 
ferent values for Rc- 

extrapolation are used in the model and shown in Fig. [4] 

i) Saturation of the correlation radius at Rc — 1 fin, 

ii) An increase of R c from SPS and RHIC to LHC. 

In the first case, namely of an energy independent Rc 
and at fixed temperature, the strangeness suppression at 
LHC will be the same as at RHIC. Consequently, different 
ratios of strange to non-strange particle yields will be 
modified only through the variation in /ib which can be 
quantified by the exp(±/^/T) Boltzmann factor. 

In the second case, of increasing Rc with y/s and at 
fixed temperature, the strangeness suppression at LHC 
will be weaker than at RHIC leading possibly to an equi- 
librated, canonical system without any additional sup- 
pression. The scenario ii) is naturally expected if Rc co- 
incides with the size of the system. In this case, the Rc 
should scale with the number of pions in the final state. 
This scenario could be verified experimentally at LHC 
by comparing the strange/non-strange particle ratios in 
p — p collisions for events with different pion multiplic- 
ities. If valid, then for sufficiently high pion multiplic- 
ity at large -y/s the strangeness production normalized to 
pion multiplicities could converge to the results expected 
in heavy-ion collisions. However, in view of the known 
data in elementary and heavy-ion collisions this is very 
unlikely scenario. An increase of Rc with yfs shown in 
Fig. [4] is expected to saturate at higher energies, or to 
be much weaker than linear. Due to lack of data the 
actual dependence of Rc = Rc{\fs) is not known. The 
LHC data are essential to understand this behavior of 
strangeness suppression and its energy dependence. 

In Fig. [5] and Table IIVI we summarize predictions of 
the statistical model for different particle ratios at LHC 



FIG. 6: Particle ratio Q/n as a function of Rc for two assumed 
temperatures, T = 160 MeV (dashed) and T = 170 MeV (full 
line) . 



TABLE IV: Particle ratios for thermal conditions expected at 
LHC. The extreme values of the cluster size Rc span the band 
of expected numerical values. The grand-canonical values in 
the last column are included for comparison, for a detailed 
discussion see Ref. [l7l |. 



Ratio 


Rc = 1 fm 


Rc = 2 fm 


grand canon. 


p/V" 


0.0970 


0.0920 


0.0914 


K+/TT + 


0.0871 


0.169 


0.180 


K~/n~ 


0.0870 


0.169 


0.179 


A/p 


0.179 


0.436 


0.473 


E-/A 


0.0397 


0.130 


0.160 


Q-/S- 


0.0358 


0.131 


0.186 


or/K~ 


2.83 • 10" 4 


4.06 • 10" 3 


7.19 • 10" 3 




2.46 • 10~ 5 


6.85 • 10~ 4 


1.29 • 10~ 3 



energy. We compare the results obtained under the GC 
description of strangeness conservation with the canon- 
ical model (c) with the parameter Rc describing the 
size where strangeness is conserved exactly in a system. 
Changing the value of Rc from one to two fermi implies 
dramatic change in the ratios involving multistrange par- 
ticles (see Fig. [5]). It is interesting to note, that even for 
Rc = 2 fm the S/V and Q/ir ratios differs from their 
GC values, whereas K/w and A/w are already well con- 
sistent with GC results. The CI/tt is particularly sensi- 
tive to changes in Rc while its temperature dependence 
is rather moderate as seen in Figures [5] and [H Conse- 
quently, the Q/tt ratio is an excellent observable to probe 
strangeness suppression and correlated strangeness pro- 
duction in p — p collisions at LHC energy. 
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V. SUMMARY 

We have analyzed, within the statistical model, re- 
cent data on particle production in p — p collisions at 
SPS and RHIC energies. The models were formulated in 
the canonical ensemble with respect to strangeness con- 
servation and extended to implement extra strangeness 
suppression either by the off-equilibrium 75 or by a 
strangeness correlation volume. We have shown, that at 
RHIC the midrapidity data can be successfully described 
by the canonical statistical model with and without any 
extra suppression effects. On the other hand, the full 
phase-space SPS data require additional strangeness sup- 
pression that is beyond the conventional canonical model. 
Extrapolating all relevant thermal parameters from SPS 
and RHIC up to LHC energy we have made predictions 
for particle yields in p — p collisions at *Js = 14 TcV. We 
have discussed the role and the influence of strangeness 



correlation volume on particle production in p-p colli- 
sions. We have argued that the fl/ir ratio is an excel- 
lent probe of strangeness correlations and/or strangeness 
suppression mechanism in p — p collisions. We have indi- 
cated that comparing strangeness production at the LHC 
in events with different pion multiplicities can provide a 
deep insights into our understanding of strangeness sup- 
pression from A — A to p — p collisions. 
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